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Abstract. Calculation of higher-order two-loop corrections is now a limiting factor 
in development of the bound state QED theory of the Lamb shift in the hydrogen 
atom and in precision determination of the Rydberg constant. Progress in the study 
,.1^ ' of light hydrogen-like ions of helium and nitrogen can be helpful to investigate these 

Mh| uncalculated terms experimentally. To do that it is necessary to develop a theory of 

such ions. We present here a theoretical calculation for low energy levels of helium 
and nitrogen ions. 
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1 Introduction 



• 1-H . The Quantum Electrodynamics (QED) theory of simple atoms like hydrogen 

►;^._^' or hydrogen-like ions provides precise predictions for different energy levels 

(~| . [0jll- Particularly, some accurate results were obtained for the Lamb shift in 

0^' the ground state of the hydrogen atom. The accuracy of the QED calcula- 

tions of the Lamb shift has been hmited by unknown higher-order two-loop 
corrections and inaccuracy of determination of the proton charge radius Q . 
^ . As far as the proton size is going to be determined very precisely from a new 

C^ ' experiment B on the Lamb shift in muonic hydrogen, the only theoretical 

uncertainty is now due to the two-loop contribution. Improvement of the the- 
ory is important to determine the Rydberg constant with high accuracy ||5|,|6| 
^~^ ' . and to test the bound state QED precisely. 

f^ ' Since the theory seems not to be able to give now any results on higher- 

^D i order two-loop corrections (a^ {Zafm and higher) we have to look for another 

c/3 ' way to estimate these terms and so the uncertainty of the hydrogen Lamb 

_^ . shift theory. An opportunity is to study the problem experimentally, mea- 

C/3 ' suring the Lamb shift in different hydrogen-like ions at not too high value of 

^^. the nuclear charge Z. Only for two such ions the Lamb shift can be available 

with a high accuracy from experiment at the present time or in the near fu- 
ture. Namely, these are helium [Q,|| and nitrogen ||^ ions. The experimental 
estimation of higher-order two-loop terms is quite of interest also because of 
^ . recent speculation on a great higher-order term [|o| (see also Refs. pi] , p^ ). 

^H ' The advantage of using Z > 1 \s determined by the scaling behaviour of 

different QED values: 
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• The scaling of the Lamb shift is Z"*; 

• The scaling of the radiative line width of excited states (e. g. of the 2p 
and 3s states) is Z^ as well; 

• The scaling of the unknown higher-order two-loop corrections to the 
Lamb shift is Z^ . 

Thus, relatively imprecise measurements with higher Z can nevertheless give 
some quite accurate data on some QED corrections. 

Our target is to develop a theory for the Lamb shift and the fine structure 
in these two atomic systems. Eventually we need to determine the 2s — 2pii2 
splitting in the helium ion (for comparison with the experiment |7l ) , difference 
of the Lamb shifts El{2s) — i?L(3s) in ^He+ (for the project ||]) and the 
2p3/2 ~ 2s interval in hydrogen-like nitrogen. The difference mentioned is 
necessary Ofl if one needs to compare the results of the Lamb shift (n — 2) 
measurement |7| and the 2s — 3s experiment. 

Since the uncertainty of the QED calculations is determined for these 
two ions (He^ and N^+) by the higher-order two- loop terms, we are going 
to reduce the other sources of uncertainty. We present results appropriate 
to provide an interpretation of the experiments mentioned as a direct study 
of the higher-order two loop corrections. The results of the ions experiments 
should afterwards be useful for the hydrogen atom. 



2 Theoretical contributions 

2.1 Definitions and notation 

The Lamb shift is defined throughout the paper as a deviation from an un- 
perturbed energy level[J 

E^^){nl,)=mn[f{nj)-l]-,^^^J^^[f{nj)~\]\ (1) 

where M and m are the mass of the nucleus and of the electron, and niB. 
stands for the reduced mass. The dimensionless Dirac energy of the electron 
in the external Coulomb field is of the form 



linj) = 



' {Zaf 



\ (,^-J-l + ^{J + lY-[ZaYy 



(2) 



The Lamb shift is mainly a QED effect, perturbed by the influence of the 
nuclear structure 

AE = Z\£;qed + AE^nci ■ (3) 

^ We use the relativistic units in which h = c= 1. 
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The QED contribution 



AE, 



QED 



= AE^ + AE 



M 



(4) 



includes one-, two- and three- loop terms calculated within the external- field 
approximation 



_ a{ZaY /mnY m 



n \ m / n"^ 



F+ {-]H+ [-] K 



(5) 



and a recoil corrections AEm, which is a sum of pure recoil and radiative 
recoil contributions depending on the nuclear mass M 



AEm = AE^ 



Rcc 



AE 



RRC 



(6) 



2.2 One-loop contributions: self energy of the electron 

Let us start with the one-loop contribution. The terms of the external-field 
contribution are usually written in the form of an expansion 



(Zc.)M^) = ^ A,,(Za)Mn^- ^ 



(7) 



«j 



The dominant contribution comes from the one-loop self energy of the elec- 
tron. The known results for some low levels are summarized below^: 



FL^iZ) 



4 m 

3 (Za)^ mji 



4^ A / ^^ 10 
-ln[ko(ns)) + — 



iZa)4.(§^-lH2))+iZaf{-ln'-J-^^ 



+ A(ji{ns) In 



(Za) 



+ Gns{Z) 



+ G2pi^2(Z) 
4 



PZ^iZ) 



-ln(ko{2p)) + —— + {Za)^ f^ln— L- 
3 ^ ^ ^'' UniR ^ ' VSO {Za)^ 



+ 62^3,2(2') 



where the state-dependent logarithmic coefficient AQi{ns) is known 

A6i(ls) = — ln(2)- — , 
^^ ' 3 ^ '^ 20' 



It is useful to keep somewhere the reduced mass itir. 



(8) 



(9) 



(10) 
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A6i(2s) = yln(2) + g, 
Aai(3,s) = -41n(3) + |ln(2) + ^, 

the Bethe logarithm is [Q 

ln(A:o(ls)) =2.984128 56... , 
ln(fco(2s)) =2.811769 89... , 
ln(fco(3s)) =2.767 663 61... , 
ln(fco(2p)) = -0.030 016 71... 

and higher-order self-energy terms Gni{Z) are nmiierically found in Sect. 2.5. 

2.3 One-loop contributions: polarization of vacuum 

The coefficient of the expansion M) for free vacuum polarization can be cal- 
culated in any order of Za in a closed analytic form p5[ . In particular, the 
result was found in Ref. |jla for the circular states {I = n— 1) with j = / + 1/2. 
For n = 1, 2 one can expand at low Za and find 

+ (Za)^ { — n- -ttH n2H , 11 

^ '^ V96 (Za)2 48 ^'^2887' ^ ^ 

Fl'' {Z) = -— [Zoif + ^^ {Zaf . (12) 

2P3/2V I 70 1024 

In the case of other states the result has been known only up to the order 
[Za)^ ||l7| , [l6| , [l3[ . We present here new results for two other states at n — 2 
and for the 3s state: 

FriZ)^-±^%iZa) 

^ ^ V 15 (Za)^ 900, 

,r^ n'^s/Stt, 1 57r , ,„, 84l7r\ , „, 

+ (Za)^ i — n- -^ + — n2)H , (13) 

^ '^ V96 (Za)2 24 ^^92167' ^ ^ 

FY/ (Z) = - — {Zaf + i^ {Zaf , (14) 

2pi/2V ) 140 3072 ^ ' ' ^ ' 

FlP{Z) = ~- + —{Za) 
3s \ ) 15 48 ^ ' 

,r, ,2 f '^ , 1 4 3 1139\ 

+ {ZaY In -rr In 

^ ^ V 15 {Zaf 15 2 1575/ 
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+ (Zaf — In 4 ln(6) H . (15) 

^ '^ V96 (Za)2 48 ^^2592/ ^ ^ 



2.4 Wichmann-Kroll contributions 

It is not enough to consider the free vacuum polarization. The relativistic 
corrections to the free vacuum polarization in Eqs. (O-H) are of the same 
order as the so-called Wichmann-Kroll term due to Coulomb effects inside 
the electronic vacuum-polarization loop. To estimate this term we fitted its 
numerical values from Ref. ||Tq], which are more accurate for some higher 
Z ~ 30, by expression ^]^ 

Fwk{Z) = {Zaf [}^-^^+ {Zaf (^An In ^ + A^o) . (16) 

We found that the contribution of Aiq and A71 terms is small enough. The 
accuracy of the calculation in Ref. |18] is not high at Z = 7 and we performed 
some fitting of higher Z data. To make a conservative estimate we find two 
pairs of coefficients which reproduce the result at Z — 30. The results are: 

An{ns) = -0.23(2), Ajo{ns) = and Fwk{7) = 0.000 139(1) ; 
An{ns) ^ 0, ^70(^5) = -0.07(1) and i^vFK (7) = 0.000 130(2) , 

where the uncertainty comes from inaccuracy of the numerical calculations 
of -^^^(30) = 0.0020 [|l^, which is estimated here as a value of a unit in the 
last digital place. Comparing the results above one can find a conservative 
estimate: FwKi7) = 0.000 134(6). The value FwKiiO) = 0.0020 is valid for 
both the Is and 2s states and we use this value for the 3s state as well. On 
this level of accuracy {6Fwk{'^0) — 0.0001) there is no shift of the 2p levels 
[Esl and we use a zero value for them. 



2.5 Fitting of one-loop self energy contributions 

We separate from the expression for the self-energy part of the one-loop 
correction (Q) the function 

G{Z) = Ago + (higher-order terms) . (17) 

Using numerical values of Aqq from Refs. ||l^,Q and ones of G{Z) from Rcfs. 
[ pip^ , we performed several types of fitting for these functions. 
We started with fitting (I) with function 

1 



Aeo + (Za) ( Aji In j^—^ + A70 ) , (18) 
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minimizing the sum 

^ / G{Z) - [Ao + {Za) (An \n{ll{Zaf) + A70)] | ' 

with respect to A70 and A-n for Is state (where Z = 1 ... 5) and to A70 
for 2s state (where Z = 5,10). In the latter case we used the fact that 
A'!i{\s) — A-ji{2s). The statistical error of data 

5G{Z)\ ' = [^„umG(Z)]' + [T:\Zaf Af^] ' . (20) 

contains uncertainty of numerical integrations in Refs. 12^,|2^ and of the fit 
in ( [l7| ) due to neglecting of higher-order terms of absolute order a{Zaf'm 

(where A\q is a result of preliminary fitting with 5G = 5numG )■ 

To estimate the additional systematic uncertainty which originates from 
the unknown term of order a(Za)^ we studied a sensitivity of the fit (I) to 
introduction of some perturbation function h{z) 

GiZ) ^ G{Z) + h{Z) , (21) 

The final uncertainty of the fit was calculated as a random sum of differences 
of the fits without function h and with h{Z) from the binomial expansion of 
the expression 

1 



3 



(^«)'(ln(^^+^) (22) 



, 2 

\2 ' 



for the Is and 2s states and 



for the 2p-states and for the difference G2S — Gns (see below). 

The logarithm ln(l/(Za)^) in the expansion (0) is a large value^ at very 
low Z but it is quite a smooth function of Z around Z = 7 (see Table p. 
Due to that, we can also use a non- logarithmic fitting function 

GiZo) + A{Z - Zo) + ^{Z - Zo)^ (24) 

with smooth behaviour at Z ^ 7. In particular, we applied the Eq. (EJ) for 
numerical data from Ref. i2^,p3| at Z — 3,4,5 with central value Zq = A (fit 
//) and Z = b, 10, 15 with Zq = 10 (fit ///). 

For the 2s state we also performed independent fits for Gis and the dif- 
ference Gis — G2S, finding G2S as their combination. Values of corresponding 

^ Note that a natural value for the constant term is about In (fco(ns)) ~ 3). 
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Table 1. Function \a{l/{Zaf) for small Z 



z 


\n{l/{Zaf) 


Z 


\n{l/{Zaf) 


1 


9.840 


7 


5.949 


2 


8.454 


8 


5.682 


3 


7.643 


9 


5.446 


4 


7.068 


10 


5.235 


5 


6.622 


15 


4.424 


6 


6.257 


20 


3.849 



fits are labeled (IV) (fits for low Z), (V) {Z = 3, 4, 5 for Is and Z = 5, 10, 15 
for tlie difference Gis - G2s) and (VI) (both fits for Z = 5J0, 15). That can 
be useful because data on the Is state is more accurate |2l|,|22| , and in case 
of difference the uncertainty is smaller (cf. Eqs . (2^ ) and (|2^)) and one of 
higher-order parameters is known (A71 = 0) [^3[|l3 ]. 

Only few data are available for 3s [gj| at Z — 10, 20, 30 . . . and we perform 
two fitting: fit / with Z = 10, 20 and fit /// with Z = 10, 20, 30 at Zq = 20. 

Different fitting functions are plotted in Figs. |^-g for n = 1, n = 2, 
Gis — G28 and G2S — G^s- The points with error bars are for numerical values 
obtained in Refs. 0,^li2l. 
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Fig. 1. Fitting of G{Z) for Is and 2s states. 



The results of all fits for Z ^ 2 are summarized in Table ^ and Z — 7 vn 
Table ||. The value for the Is state at Z = 7 is in agreement with the less 
accurate one obtained in Ref. ESll. 



8 Vladimir G. Ivanov and Savely G. Karshenboim 



G(1s)-G(2s) 




0,05^ 




G(2s)-G(3s) 










(III) 






^^^^ 




0,10- 
0,15- 









3 6 9 12 15 18 21 24 27 30 



Fig. 2. Fitting of Gi,(Z) - G^,{Z) and G-2.s{Z) - Gi.lyZ). 
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Fig. 3. Fitting of G(^) for the 2pi/2 and 2p3/2 states. 



2.6 Two-loop contributions 



Two-loop corrections have not yet been calculated exactly and only a few of 
the terms have been known up-to-date pfl.Bl[: 



HnsiZ) 



w ~ i"' + r' ^"(^) - ^(^^ ' ^^-^^^(^^ (^") 



(Za)^ 



In' 



1 



27 (Za)2 



+ G^U^) 



H2p,,^{Z) = - -ai + {Zaf 



9 



-1^ (z^ + ^^--(^) 



^^'(z!^ + ^^-^(^) 



(25) 
(26) 
(27) 
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Table 2. Values G(2) for different types of fitting 



Value 


(I) 


(IV) 


Gls 


-29.772(7) 




G2s 


-30.7(2) 


-30.64(3) 


Gls — G2s 


0.87(3) 




G2s — G3.S 


-0.12(9) 




<-'2pi/2 


-0.95(3) 




^2^3/2 


-0.48(3) 





Table 3. Values G(7) for different types of fitting 



Value 



(J) 



(W 



(III) 



(IV) 



(V) 



(VI) 



Gls 


-27.6(1) 


-27.68(1) 


-27.677(9) 


G2s 


-28.5(1) 




-28.47(5) 


Gls — G2s 


0.77(3) 




0.79(5) 


G2pl/2 


-0.84(3) 




-0.85(5) 


G2p3/2 


-0.41(3) 




-0.41(3) 



-28.4(1) -28.47(5) -28.45(3) 



where for I ^ the leading contribution is due to the anomalous magnetic 
moment of the electrory 



197 7r2 7r2 3 

«4 = ^ + ^-yln(2) + -C(3) 



-0.328 478 966. 



(28) 



The higher-order terms denoted as G^^{Z) are not known. We expect |Q 
that the magnitude of G^^ does not exceeded by half the value of the leading 
logarithmic term (In^(Za) for the ns states and \v?{Za) for the 2p states) 
at Z — 1 and use the estimation for G^' {1} for any Z. The purpose of this 
calculation is to eliminate any other sources of theoretical uncertainty and to 
study a value of 6*2^ by a comparison with experiment. 

The leading state-dependent term for the s-state is also known ||3^ , |l3| , 
and in particular one can find 



H,s{Z) - H2s{Z) = (Za)^ 
H2s{Z)-HUZ) = {Zaf 



In 



In^ 



2 1 ^16,_^,^, 7^ , ^ji 



1 /16, 3 91 > ^Ilrry. 

(Z^U^^2-8t)+^-(^) 



,(29) 
,(30) 



We follow a notation {g^ — 2)qed = 2 • ^ a2n (o/tt) . See Ref. M for detail. 
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and unknown higher-order contributions G^^ are estimated in the same way 
as for the 2p states. 



2.7 Three-loop contributions 

The three-loop contributions for ns state at Z — was eventually obtained 
in Ref. §^ 

X„,(0) = 0.4174..., (31) 

while for the p-states the corrections 



■f*^2pi/2(0) =" ^3^6. 



K. 



2p3/2 



(0) 



0-6 



comes from the g — 2oi electron: 

ag = 1.181241... 



(32) 
(33) 

(34) 



2.8 Pure recoil corrections 

The pure recoil correction 



■K n'' M 



where is known analytically with sufficient accuracy |26 
Rns{Z) = [Za 



2, 1 8, ., ... 187 



(Za)^^(41n(2)-I) , 



R2j>{Z) = {Zc 



--ln(fco(2p))-- 



{Zaf - 



(35) 



(36) 

(37) 



Some numerical results for the states with n — 1,2 are also available |^ . 
Second-order recoil corrections {{Za)'^{m^/AP)) are known |g8[ but their 
contribution is negligible. 



2.9 Radiative-recoil corrections 

The radiative-recoil corrections are known only in the leading order 

AEjiRcinl) = ^ '- -^ — -1.36449 Sm . 



(38) 
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2.10 Finite-nuclear-size correction 

The correction for s-states due to the finite size of the nucleus is of the form 

/ii?Nuci(nO - l^-^^P^ {mR^ySio (l + (^")'li^^^;^) ' (39) 

where Rn is the rnis nuclear charge distribution radius. The leading term in 
Eq. (p^) vanishes for the p states because their non-relativistic wave function 
is vanishes at the origin itself. However, the small component of the Dirac 
wave function contains a factor of ap/2m and the Dirac wave function is not 
equal to zero at the origin. In particular, one can easily find 

2 3 {Za)^m^ 



32 



and 



('^2p3/.(0))' = 0, (40) 

Z\^Nuci(2pi/2) = — (Za)^ m (m T^n) , 

Zi^Nucl(2p3/2) = . (41) 

The error of the logarithmic term can be estimated as 1/2 Z\_Enuc1 (Za)'^ ln(l/a). 
The nuclear radii and corrections are presented in Table 0. 



Table 4. Nuclear square charge radius 



Atom 


Rn 


Z\£nuci(2s) 


Ref. 




[fm] 


[MHz] 




^He 


1.674(12) 


8.80(12) 


§1 


"N 


2.560(11) 


3145(33) 


fei 


15n 


2.612(9) 


3274(30) 


© 



3 Summary 

All contributions and final values of the Lamb shifts in the lowest states of 
hydrog-like ions '*He+, ^''N^"'' and i^N^+ are presented in Tables H, and M, 
respectively. The final results for intervals which can be measured are listed 
in Table @. The results in Table involve three main sources of uncertainty 
and we split the uncertainty there and in auxiliary Tables &fl, respectively: 
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• The higher-order two-loop corrections; 

• Nuclear structure corrections (They require further study); 

• Other theoretical QED uncertainties beyond the higher-oder two-loop 
effects. 

The theory (with unknown higher-order two-loop effects excluded) is found to 
be accurate enough and we hope the study of helium and nitrogen hydrogen- 
like ions is a promising way to study in detail the two-loop contributions 
experimentally. 



Table 5. Different contributions to the Lamb shift in hydrogen-like helium *He (in 
MHz) 



Contribution 



2s 



2pi/2 



2p3/2 



SE 


14 251.51(2) 


-204.747(10) 


201.503(10) 


VP 


-426.78 


-0.022 


-0.005 


WK 


0.02 








One-loop 


13 824.75(2) 


-204.769(10) 


201.498(10) 


Two-loop 


0.70(11) 


0.420(2) 


-0.201(2) 


Three-loop 


0.004 


0.003 


0.002 


Pure recoil 


2.53 


-0.131 


-0.131 


Radiative recoil 


-0.01 









QED 
Nuclear size 



13 827.98(11)(2) 
8.80(13) 



-204.484(2)(10) 




201.168(2)(10) 




Total 



13 836.77(11)(13)(2) -204.484(2)(10) 201.168(2)(10) 
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